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1.  Introduction 


Magnetically  soft  (high  penneability)  materials  are  often  used  as  magnetic  field  concentrators  in 
magnetic  sensor  systems.  These  magnetic  flux  concentrators  improve  the  sensitivity  of  the 
associated  system  by  increasing  the  magnetic  flux  at  the  position  of  the  sensor.  Flux 
concentrators  have  been  used  with  spin  valves  (1 ),  hall  sensors  (2),  atomic  magnetometers  (3), 
superconducting  quantum  interference  devices  (SQUIDs)  (4),  resonating  magnetic  sensors  (5), 
and  magnetic  tunnel  junctions  ( 6 ).  The  concentration  effect  is  a  function  of  (1)  the  relative 
penneability  (jUr)  of  the  material  used  to  make  the  flux  concentrators,  (2)  the  area  ratio  between 
the  ends  of  a  concentrator  (7),  and  (3)  the  size  of  the  air  gap  between  adjacent  flux  concentrators 
(5).  It  has  been  shown  that,  with  judicious  choices  of  flux  concentrator  geometry  and  air  gap 
size,  magnetic  field  amplification  can  exceed  the  theoretical  value  fixed  by  jur  of  the  material 
used  (7). 

Magnetoresistive  (MR)  devices,  particularly  magnetic  tunnel  junctions  (MTJs),  are  excellent 
candidates  for  low-cost,  low-power  magnetic  sensors.  While  MTJs  can  provide  excellent 
detectivity,  taking  full  advantage  of  these  values  at  low  frequencies  (f  <  100  Hz)  requires  finding 
a  method  of  mitigating  1// noise  (9, 10).  We  have  reported  elsewhere  (11)  on  the  successful 
fabrication  and  validation  of  a  concept  for  mitigating  the  effect  of  1  If  noise  by  increasing  the 
operating  frequency  of  the  magnetic  sensor  to  a  region  where  1  If  noise  is  much  smaller.  The  idea 
is  to  place  the  magnetic  sensor  between  flux  concentrators  that  have  been  deposited  on 
microelectromechanical  system  (MEMS)  flaps.  An  essential  part  of  designing  the  device  was 
understanding  the  behavior  of  the  magnetic  field  lines  both  for  enhancing  the  field  at  the  position 
of  the  sensor  as  well  as  maximizing  modulation  of  the  field.  These  needs  motivated  our  work  on 
macro-magnetic  modeling. 


2.  Concept  of  the  MEMS  Flux  Concentrator 


The  key  concept  of  the  MEMs  flux  concentrator  is  that  the  flux  concentrators  are  set  into 
oscillatory  motion.  Once  set  into  periodic  motion  at  frequencies  of  around  10  kHz  or  more,  the 
enhancement  of  the  magnetic  field  at  the  position  of  the  sensor  is  modulated.  Due  to  this  high 
frequency  modulation,  the  sensor  operates  at  a  frequency  where  the  1  If  noise  is  much  smaller. 
Sidebands  to  the  drive  signal  and  its  first  harmonic  appear  due  to  this  modulation.  There  are  at 
least  two  major  matters  that  have  to  be  addressed  before  the  potential  advantages  of  this  device 
can  be  realized.  First,  it  is  necessary  that  the  sensor  element  is  responsible  for  most  of  the  1  If 
noise  and  not  some  other  part  of  the  sensor  system.  The  second  issue,  the  design  and  fabrication 
of  the  device,  is  discussed  in  the  sections  concerning  specific  designs. 
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As  described  here,  the  MEMS  flux  concentration  will  greatly  reduce  the  1  If  noise  of  the  magnetic 
sensors.  It  does  nothing  to  reduce  the  1  If  noise  of  the  flux  concentrators  themselves.  Thus,  if  the 
concept  is  going  to  be  useful,  it  is  necessary  that  the  1  If  noise  of  the  flux  concentrators  be  much 
less  than  the  1  If  noise  of  the  magnetic  sensors.  This  requirement  motivated  a  series  of 
experiments  to  see  if  this  requirement  was  fulfilled.  As  we  reported  elsewhere  (12),  the  noise 
voltage  of  spin  valves  with  and  without  adjacent,  stationary  flux  concentrators  are 
indistinguishable.  This  result  implies  that  the  1  If  of  the  flux  concentrators  is  much  less  than  the 
1  If  noise  of  the  spin  valves. 


3.  Magnetic  Modeling  Details  and  Parameters 


We  perfonned  magnetic  modeling  using  a  commercial  finite  element  code  program  called 
Maxwell  3D,  from  ANSYS.  Maxwell  3D  is  capable  of  analyzing  AC  magnetic,  DC  magnetic, 
and  electrostatic  field  problems.  The  three-dimensional  (3-D)  DC  magnetic  portion  of  the 
software  computes  static  magnetic  fields,  where  the  source  originates  from  a  DC  current  or 
voltage,  pennanent  magnets,  or  externally  applied  magnetic  fields.  It  can  directly  compute  the 
magnetic  field  (H)  and  current  distribution  (J),  and  derive  the  magnetic  flux  density  (B)  from  the 
H  field.  In  addition,  it  can  automatically  calculate  force,  torque,  inductances,  and  saturation  in 
devices  containing  linear,  nonlinear,  and  anisotropic  materials.  The  post-processor  portion  of  the 
software  can  provide  plots  of  flux  lines,  B  and  H  fields,  energy  densities,  and  saturation.  The 
modeling  process  consists  of  drawing  the  objects  of  interest,  assigning  properties  (coercivity, 
penneability,  etc.)  to  the  objects,  assigning  boundaries  or  sources,  seeding  the  objects  and 
creating  a  mesh,  and  then  the  processing  of  the  now  defined  problem. 

Each  model  investigated  involved  the  same  fundamental  sequence  of  steps.  The  first  step  in  the 
analysis  of  the  flux  concentrator  was  to  draw  the  model.  Drawing  the  model  consists  of  drawing 
3-D  objects  and  either  joining  them  together  or  subtracting  them  from  each  other.  This  allows 
one  to  create  complex  objects.  A  sufficiently  large  region  around  the  flux  concentrators  and 
sensing  region  was  defined  as  a  background.  The  flux  concentrators  are  drawn  as  solid  pieces 
and  assigned  the  material  properties  of  permalloy  (nickel-iron  [NiFe])  with  a  penneability  of 
5,000  as  that  is  a  value  that  is  readily  achieved  in  thin  films  of  the  material.  The  material 
properties  assigned  to  this  background  are  those  of  a  vacuum,  with  a  relative  permeability  of  1 . 
The  modeling  is  macro-magnetic  in  nature  as  it  does  not  take  into  account  domain  structure  but 
does  incorporate  demagnetization  factors.  The  initial  mesh  is  created  by  the  program  but  one  can 
create  regions  in  which  the  initial  mesh  is  denser  so  as  to  force  more  tetrahedrons  into  regions 
where  one  has  a  greater  interest  in  the  solutions  without  significantly  increasing  solution  time. 
Mesh  refinement  is  also  handled  by  the  program  as  part  of  an  iterative  process  in  which  energy 
error  and  percentage  decreases  to  a  predetermined  figure. 


2 


In  order  to  determine  the  accuracy  of  the  modeling  results,  measurements  were  made  using  spin 
valves  and  fixed  flux  concentrators.  Two  different  sizes  of  trapezoidal-shaped  flux  concentrators 
were  fabricated  and  measurements  of  the  magnetoresistance  were  made  for  (1)  only  spin  valves, 
(2)  a  spin  valve  between  a  pair  of  “normal”  flux  concentrators,  and  (3)  a  spin  valve  between  a 
pair  of  larger  flux  concentrators.  The  smaller  pair  of  flux  concentrators  had  a  height  of  104  pm, 
a  short  side  length  of  80  pm,  and  a  long  side  length  of  150  pm.  The  larger  pair  of  flux 
concentrators  was  80  pm  long  on  the  short  side  and  300  pm  long  on  the  long  side,  and  had  a 
height  of  3 15  pm.  The  air  gap  between  both  types  of  flux  concentrators  was  held  constant  at 
52  pm.  Figure  1  shows  the  resistance  of  the  spin  valves  with  and  without  the  flux  concentrators. 
The  amount  that  the  flux  concentrators  enhanced  the  field  was  compared  to  calculated  values 
from  magnetic  modeling  and  the  results  are  shown  in  table  1 ,  showing  that  the  measured  field 
enhancement  provided  by  stationary  flux  concentrators  agrees  to  within  4%  of  the  values 
detennined  from  the  model  results.  These  results  are  particularly  good  considering  the  flux 
concentrators  were  modeled  as  solid  pieces  of  permalloy,  while  the  actual  samples  consisted  of 
two  repeats  of  40  A  chromium  (Cr)/1500  A  pennalloy. 


Figure  1.  (a)  Magnetoresistance  of  a  spin  valve  with  small  and  larger  flux  concentrators  and  (b)  an  enlarged  view  of 
the  data  near  zero  applied  field. 
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Table  1.  Enhancements  provided  by  flux  concentrators. 


Ave.  Linear 
Region  Slope 
(mV/Oe)* 

Measured 

Enhancement 

Calculated 

Enhancement 

No  Flux 
Concentrator 

0.21133 

- 

- 

Small  Flux 
Concentrator 

0.41826 

1.98 

1.93 

Large  Flux 
Concentrator 

0.69586 

3.29 

3.36 

*Voltage  across  spin  valve  was  set  to  ~400  mV  in  all  cases. 


We  are  interested  in  the  magnetic  field  strength  and  flux  line  behavior  at  the  position  the  sensor 
would  occupy.  Two  main  quantities  of  interest  to  us  are  the  enhancement  factor  and  the 
percentage  modulation.  The  enhancement  factor  EH  is  defined  as  Hs/Happi,  where  Hs  and  Happi 
denote  the  magnetic  field  strength  at  the  position  the  sensor  would  occupy  and  the  magnetic  field 
strength  far  from  the  concentrators,  respectively.  For  the  rotating  flux  concentrator  design  and 
the  out  of  phase  design,  we  define  the  percentage  of  modulation  of  the  field  at  the  position  of  the 
sensor  as 


[(HS2-Hsi)/Hsr]*  100,  (1) 

where  Hs2  and  Hsi  are  the  enhancement  factors  at  the  two  end  points  of  motion  and  Hsr  is  the 
enhancement  factor  when  the  flux  concentrators  are  at  rest  in  their  neutral  position.  Detennining 
the  modulation  for  the  in-phase  design  is  discussed  in  the  corresponding  section.  The  next 
sections  discuss  the  three  primary  designs  we  have  modeled. 

3.1  Pivoting  Flux  Concentrators 

The  first  design  of  the  MEMS  flux  concentrator  used  electrostatic  plates  to  drive  the  motion  of 
the  concentrators.  The  concentrators  perform  oscillatory  rotation  about  the  torsional  suspensions 
(figure  2),  thus  modulating  the  field  at  the  position  of  the  sensor.  Using  the  geometry  of  the 
“normal”  concentrators  flaps,  described  previously,  and  assuming  both  a  3 -pm  gap  and  0.5-pm 
thickness  of  permalloy,  we  ran  models  to  detennine  how  the  enhancement  of  the  field  at  the 
center  of  the  air  gap  would  vary  with  penneability.  One  sees  in  figure  3  that  for  permeabilities 
near  the  values  expected  for  permalloy  (around  5000)  the  enhancement  is  a  rather  weak  function 
of  the  permeability.  Thus,  it  is  not  necessary  to  maximize  the  permeability.  Holding  the 
penneability  constant  we  then  looked  at  how  rotation  of  the  flaps  alters  the  enhancement  for 
various  air  gap  sizes.  For  a  given  gap  size,  the  model  was  run  for  each  degree  of  tilt  angle  from 
0°  to  6°.  The  results  are  shown  in  figure  4  and  one  sees  that  a  large  modulation  of  the  field  can 
be  achieved  at  the  position  of  the  sensor  though  the  modulation  decreases  rapidly  as  the  distance 
between  the  flaps  is  increased.  Unfortunately,  the  processing  necessary  to  liberate  the  flaps  for 
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rotation  consistently  destroyed  the  magnetic  sensor  stack,  regardless  of  what  sort  of  protective 
coating  was  used.  This  necessitated  a  change  in  design. 


Figure  2.  Schematic  of  the  rotating  MEMS  flux  concentrator.  All  dimensions  are  in  microns. 
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Figure  3.  Magnetic  field  enhancement  at  the  center  of  the  air  gap 
when  the  flaps  are  in  the  plane  of  the  sensor  position. 
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Figure  4.  Enhancement  of  the  magnetic  field  at  the  position  of  the 
magnetic  sensor  as  a  function  of  tilt  angle  of  the  MEMS 
flaps  for  various  air  gap  sizes. 


3.2  Comb  Drive,  Out-of-phase  Flux  Concentrators 

To  avoid  damaging  the  magnetic  sensors  during  the  release  step  of  fabricating  the  MEMS 
structure  by  exposing  them  to  hydrofluoric  acid  (HF),  two  chips  were  used:  one  with  the  MEMS 
structure  and  the  flux  concentrators  and  the  other  with  the  magnetic  sensors.  The  two  chips  were 
combined  with  flip  chip  bonding.  The  sensor  chip  was  separated  from  the  MEMS  chip  by  about 
4  pm.  The  MEMS  spring  widths  were  varied  between  4  and  10  pm.  The  spin  valves  on  the 
underneath  side  of  the  sensor  chip  are  12  pm  wide.  The  two  MEMS  flaps  on  each  side  of  the 
spin  valve  were  connected  by  silicon  springs  so  that  the  desired  in  plane,  out-of-phase  motion  of 
the  flaps  was  a  nonnal  mode.  The  MEMS  flaps  are  driven  to  move  at  this  nonnal  mode  by 
electrostatic  comb  drives.  It  was  determined  that  for  proper  release  of  the  concentrators,  etch 
holes  would  need  to  be  added  to  the  flaps.  Further  details  on  this  design  have  been  previously 
published  (11,  12). 

Three  different  sets  of  magnetic  modeling  were  perfonned  for  this  version  of  the  MEMS  flux 
concentrator:  (1)  solid  flaps,  (2)  flaps  with  slit  shaped  release  holes,  and  (3)  flaps  with 
squareshaped  release  holes.  In  all  of  these  models,  the  field  values  were  determined  at  the  new 
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position  of  the  sensor,  4  pm  above  the  plane  of  the  flaps.  The  flaps  were  modeled  as  solid 
pennalloy  and  all  nonmagnetic  structural  material  was  ignored.  The  magnetic  material  of  the 
spin  valve  was  also  ignored  as  the  thicknesses  of  the  various  layers  are  small  enough  to  have 
only  a  small  influence  on  the  surrounding  flux  environment.  Enhancement  curves  for  each  case 
are  created  from  the  results  of  running  models  for  various  air  gap  sizes  (figure  5).  The  out-of¬ 
phase  motion  causes  the  air  gap  to  decrease  (increase)  as  the  flaps  move  towards  (away)  from 
each  other.  The  enhancement  curve  is  used  to  determine  the  percentage  field  modulation  at  the 
position  of  the  sensor  based  on  how  far  the  flaps  move.  Equation  1  and  the  results  of  the  models 
for  the  square  etch  holes  were  used  to  create  a  modulation  curve  (figure  6).  From  the  modulation 
curve,  we  were  able  to  detennine  that  achieving  50%  modulation  of  the  magnetic  field  at  the 
position  of  the  sensor  would  require  38  pm  of  total  travel  or  ±19  pm  of  motion  of  the  flux 
concentrators. 


Figure  5.  Enhancement  curves  generated  from  models  of  permalloy  flux 
concentrator  pairs  with  and  without  release  holes. 
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Figure  6.  Modulation  as  a  function  of  the  motion  of  the  flux  concentrator  flaps. 


In  an  effort  to  increase  not  only  the  enhancement  of  the  field  at  the  position  of  the  sensor  but  also 
increase  the  modulation  beyond  a  factor  of  2  without  needing  19  pm  of  motion,  we  modeled  the 
addition  of  a  second  pair  of  larger,  stationary  flux  concentrators.  These  stationary  flaps  were  to 
be  on  the  wafer  containing  the  sensor  and  were  drawn  to  have  a  height  of  2000  pm,  a  short  side 
length  of  125  pm,  and  a  long  side  length  of  150  pm.  The  thickness  of  the  permalloy  for  these 
concentrators  was  set  at  0.5  pm  and  they  had  a  50%  overlap  of  the  moving  concentrators.  As 
shown  in  figure  7,  the  field  enhancement  at  the  position  of  the  sensor  was  indeed  increased  but 
the  modulation  was  still  roughly  a  factor  of  2. 
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Figure  7.  Enhancement  curves  demonstrating  effect  of  the  addition  of  a  second  pair 
of  large,  stationary  flux  concentrators. 


3.3  Comb  Drive,  In-phase  Flux  Concentrators 

Since  the  magnetic  modeling  indicated  that  the  addition  of  stationary  flux  concentrators  had  little 
effect  on  the  modulation,  another  fundamental  change  was  made  to  the  overall  design.  The  main 
changes  were  (1)  the  air  gap  was  decreased  to  1 1  pm,  (2)  the  size  of  the  moving  concentrators 
was  increased,  (3)  the  number  and  type  of  etch  holes  was  increased,  and  (4)  the  moving 
concentrators  were  joined  together  by  small  struts  to  either  side  of  the  air  gap  (figure  8a).  These 
changes  allowed  us  to  change  the  nonnal  operation  mode  to  an  in-phase  motion  that  would 
alternately  enhance  the  field  at  the  position  of  the  sensor  or  shunt  the  field  away  from  the  sensor, 
thus  increasing  the  percentage  modulation.  The  flaps  now  had  a  height  of  1 10  pm,  a  short  side 
length  of  80  pm,  and  a  long  side  length  of  308  pm.  Smaller,  diamond  shaped  etch  holes  were 
placed  near  the  air  gap  to  facilitate  full  release  of  the  flaps.  The  flaps  were  modeled  as  solid, 

0.25  pm  thick  permalloy  and  all  nonmagnetic  structural  material  was  ignored  (figure  8b).  In 
order  to  simplify  shunting,  the  position  of  the  sensor  remains  roughly  4  pm  above  the  moving 
flaps  but  is  shifted  laterally  so  as  to  be  16  pm  away  from  the  air  gap  edge  of  one  of  the  flaps. 
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The  enhancement  curve  was  generated  from  the  magnetic  models  and  the  corresponding 
modulation  curve  for  the  new  in-phase  motion  is  shown  in  figure  6.  Equation  1  is  still  used  to 
detennine  percentage  modulation  but  now  both  Hsl  and  Hsr  are  reduced  due  to  shunting  and  the 
sensors  location  relative  to  the  flaps.  We  see  that  now  only  12  pm  of  motion  would  result  in 
nearly  350%  modulation  of  the  field  at  the  position  of  the  sensor. 


Figure  8.  (a)  Optical  microscope  image  of  the  in-phase  flux  concentrator  design  and  (b)  the  concentrators 
as  drawn  for  modeling.  The  darker  trapazoids  are  part  of  the  stationary  concentrators,  the 
remainder  being  out  of  view. 

The  new  location  of  the  sensor  necessitated  us  running  additional  models  to  detennine  if  the  flux 
environment  around  the  sensor  would  create  problems.  Specifically,  we  were  concerned  with  the 
flux  line  behavior  around  the  release  holes.  As  the  flux  lines  will  both  flow  around  the  holes  and 
arc  out  above  and  below  them  (figure  9),  we  wanted  to  see  if  5  pm  above  release  holes  would 
place  the  sensor  in  a  noisy,  turbulent  flux  environment.  We  see  in  figure  10  that  the  flux  lines 
are  indeed  flowing  around  the  release  holes.  The  rectangular  regions  to  either  side  of  the  flaps 
indicate  a  decrease  in  field  lines  because  this  plot  shows  field  in  the  plane  of  the  moving  flaps 
and  those  regions  are  where  the  field  is  largely  in  the  larger,  stationary  flaps  that  overlap  the 
moving  flaps  (figure  8b).  We  were  able  to  detennine  that  the  flux  environment  above  the 
moving  flaps  is  relatively  smooth  until  one  is  1  pm  away  from  the  flap  surface  (figure  11).  Only 
then  do  the  models  indicate  variations  in  the  magnetic  field  strength  due  to  the  etch  holes. 
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Figure  9.  (a)  Top  view  and  (b)  side  view  of  magnetic  flux 

lines  around  concentrators.  Color  from  blue  to  red 
denotes  low  to  high  flux  density,  respectively. 


Figure  10.  Flow  of  magnetic  flux  lines  around  etch  holes 
in  concentrators. 
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Figure  11.  Variations  in  the  magnetic  field  strength  at  1  micron  above  concentrators  due  to  etch  holes. 


4.  Discussion  and  Conclusions 


We  have  described  the  magnetic  modeling  that  has  facilitated  the  development  of  the  MEMS 
flux  concentrator,  a  device  that  uses  magnetic  flux  concentrators  deposited  on  MEMS  structures 
that  modulate  low  frequency  signals  at  the  position  of  a  magnetic  sensor.  We  presented 
magnetic  modeling  results  of  various  designs,  including  our  current  design,  focusing  on  key 
design  elements.  A  comparison  of  measured  magnetic  field  enhancement  values  to  values 
calculated  from  magnetic  modeling  was  presented  and  shown  to  agree  to  within  4%.  Working 
devices  have  been  packaged  for  testing.  In  figure  12,  we  see  four  of  our  devices  ready  for 
testing.  The  four  etched  windows  allow  us  to  verify  motion  of  the  devices  as  well  as  determine 
the  resonant  frequency  of  each  device.  As  described  in  section  3.2,  these  devices  are  flip  chip 
bonded  with  the  MEMS  devices  on  the  bottom  wafer  piece  and  the  magnetic  sensing  elements 
are  on  the  top  piece,  centered  over  the  gap  of  each  pair  of  flux  concentrators.  Using  this  packing 
we  have  been  able  to  confirm  the  concept  of  the  MEMs  flux  concentrator  (11). 
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Figure  12.  Four  MEMs  flux  concentrator  devices  packaged  for  testing. 

The  use  of  macro-magnetic  modeling  has  allowed  us  to  gain  insight  in  how  design  features 
impact  magnetic  flux  lines,  and  achieve  both  high  field  enhancement  and  a  large  percentage  field 
modulation.  Work  is  continuing  on  the  in-phase  design  and  it  is  anticipated  that  vacuum 
packaging  of  the  devices  and  the  associated  increase  in  Q,  the  resonant  frequency  divided  by  the 
frequency  halfwidth  of  the  resonant  frequency,  combined  with  the  use  of  very  sensitive  magnetic 
tunnel  junctions,  will  result  in  a  detectivity  of  a  few  pT/Hz  at  1  Hz. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


3-D 

three-dimensional 

ARL 

U.S.  Army  Research  Laboratory 

Cr 

chromium 

HF 

hydrofluoric  acid 

MEMS 

microelectromechanical  system 

MgO 

magnesium  oxide 

MTJ 

magnetic  tunnel  junctions 

NiFE 

nickel-iron 

SQUIDs 

superconducting  quantum  interference  device 
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